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Abstract 

We present analytic expressions for the vacuum polarization effects due to a light 
fermion with finite mass in the binding energy and in the wave function at the origin 
of QED and (weak coupling) QCD non-relativistic bound states. Applications to 
exotic atoms, T(ls) and ti production near threshold are briefly discussed. 
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1 Introduction 



It is well known that the vacuum polarization effects due to light fermions produce leading 
corrections to the Coulomb-like behavior of non-relativistic bound states both in QED [1] 
and in QCD [2]. If a light particle has a mass {mi) of the order of the inverse Bohr radius, 
it can neither be approximated by a massless particle nor by a very heavy one. Hence any 
observable related to the bound state depends non-trivially on mi. For a given bound state 
the mi dependence in physical observables is usually calculated numerically [1,3-8] and only 
a few analytical results are available [9-12]|3]. In this letter we present further analytical 
formulas with the exact light fermion mass dependence for the leading corrections to the 
energy shift for arbitrary quantum numbers (n, /) and to the wave function at the origin for 
the ground state. 

These formulas may be useful for quite a few physical systems of current interest. Any QED 
bound state built out of particles heavier than the electron may require them. For instance 
di-muonium, muonic hydrogen, pionium, pionic hydrogen and other simple hadronic atoms 
where the electron mass (mg) is such that me ~ fia/n, fi and n being the reduced mass 
and the principle quantum number respectively. It is worth mentioning that simple hadronic 
atoms are experiencing a renewed interest because they may allow to extract important 
information on the QCD scattering lengths for several isospin channels [13]. In particular 
the measurement of the decay width of pionium [14] at the 10% level will allow to extract a 
combination of scattering lengths with sufficient accuracy as to discern between the large and 
the small quark condensate scenario of QCD, namely between Chiral Perturbation Theory 
[15] and Generalized Chiral Perturbation Theory [16]. In fact, the leading corrections to the 
pionium decay width have been recently calculated in a systematic way using non-relativistic 
effective field theory techniques [9,17] (see [18], [5,19] and [20] for earlier relativistic, non- 
relativistic, and quantum mechanical calculations respectively). Fully analytic results have 
been obtained except for the contribution of the electron mass to the vacuum polarization 
where only a numerical result is available [5]. We shall fill this gap here and present the 
only remaining piece to have the full leading corrections to the pionium decay width in an 
analytic form. 

For QCD non-relativistic bound states, T(ls) seems to be the only one amenable to a weak 
coupling analysis [21]. Using the results of [9], analytic expressions for the binding energy 
shift due to the finite charm mass have been recently presented in [22][^ We present here 
analytic results for these effects in the wave function at the origin. 

An important non-relativistic weak coupling QCD system for the Next Linear Collider 
physics is the top quark- ant iquark pair near threshold. The production cross-section has 
already been calculated at NNLO [24]. However, the calculations are done assuming the 
mass of the bottom and charm quarks zero. Our results may allow to include the leading 

^ We learnt about the three last references after completing our calculations. 

^ The extraction of the bottom MS mass from the T(ls) mass also requires the bottom pole mass 

dependence on rric [23]. 
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effects of these masses. 



The fact that our results can be apphed to such a variety of a priory quite different physical 
systems can be easily understood in terms of modern effective field theories. Non-relativistic 
bound states have at least three dynamical scales: the hard scale (mass of the particles 
forming the bound state), the soft scale (typical relative momentum in the bound state) 
and the ultrasoft scale (typical binding energy in the bound state). Upon integrating out 
the hard scale local non-relativistic effective theories arise. These are NRQED for QED, 
NRQCD for QCD [25], and NR^L for the Chiral Lagrangian[^ [26]. Upon integrating out 
the soft scale (see [28,29] for the precise statements) effective theories which are local in 
time but non-local in space arise. The non-local terms in space are nothing but the usual 
quantum-mechanical potentials and only ultrasoft degrees of freedom^ are left dynamical. 
The corresponding non-relativistic effective theories have been named pNRQED, pNRQCD 
and pNR^L for QED, QCD and the Chiral Lagrangian respectively [26,30], (the "p" stands 
for potential). Since the leading (mass independent) coupling of the photon field to the 
non-relativistic charged particles as well as the one of the gluon field to the non-relativistic 
quarks in the NR theories is universal, it produces the same potential in the pNR theories, 
and hence they all can be discussed at once. If there is a light (relativistic) charged particle 
in QED or a light (relativistic) quark in QCD whose mass is of the order of the soft scale, it 
must be integrated out keeping the mass dependence exact, which produces a light fermion 
mass dependent correction to the static potential. 




Fig.l: Matching between the non-relativistic theory and the potential one. 



When matching the NR theories to the pNR theories only the diagram of Fig. 1 gives rise 
to a potential which contributes to the leading effect. For QED (on-shell scheme) it reads, 

Kpc(|x|) = — — / dv^ vCT (1) 

' ' ^ ' 



and for QCD (MS): 

^ We apologize to the authors of ref. [5] for slightly modifying their previously given name, 
namely Non Relativistic Chiral Perturbation Theory. The reason is that the calculation in the 
non-relativistic regime cannot be organized according to the chiral counting any longer [9]. 

In the language of the threshold expansions in QCD [27] these correspond to ultrasoft gluons 
and potential quarks. 



2 



mt 



TT |x| J [1 —V^) 6 \ 



.0 

A^^ - 1 4 1 

If N f is the number of flavors lighter that m;, the OLgiv) above runs with + 1 flavors. 
Notice that the difference between the QED and the QCD case is, apart from the trivial 
color factors «/|x| C^a^/lxl and aj-n Tfols/tt-i a term which can be absorbed in a 
redefinition of the Coulomb potential [31]. Hence for the actual calculation we shall only 
deal with (1) and use these facts to extend our results to the QCD realm. 



2 Energy Shift 



For the energy shift we obtain := 



5Er,m = -^£^n - Y< + (^(2^ + l) + {n + l){n-l- 1)) e 
-Trn Q(n + l)(2n + 1) + (n + /)(n - Z - 1)) C^- 



(2n-l)! \ ^ j ^ 2Z + 1 + A; 



2n- 



, arccos \ if ^ > 1, 
1 if e = 1, 



if e < 1. 



where — —iio?l2v? is the Coulomb energy. For ^ large, namely mi » /la/n, it reduces 



to 



2a E„ 



[n + iy. 



(2/)! 



Stt [ (n - / - 1)1(2/ + 1)! {21 + 1)!! 
^ (2/ + 2) (2/ + 2)(2/ + 4) 

, (2/ +3) (2/ +3) (2/ +5)"^ Uyyj^ 



(4) 
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whereas for ^ small, namely mi « ^a/n, we obtain 



^ 0) 



—Tin 




(5) 



where we have used (13). The key steps to obtain (3) are given in the Appendix B. We have 
done the following checks. For the IS state (3) reduces to the formula (5.3) of ref. [9]. The 
energy shifts for the IS, 2S, 2P, 3S, 3P and 3D states agree with the early analytical formulas 
of ref. [10]. For ^ large, we reproduce the well-known positronium like limit for I = (to be 
precise we agree with the correction to the energy obtained using formula (2.8) of [28]). We 
also agree for Z = 1 with formula (32) of ref. [4]. For C, small, we can compare with known 
results for massless quarks in QCD. For arbitrary n and / we agree with formula (13) of 
ref. [31]. For / = n — 1 we agree with 0(^°) and 0{C,^) of formula (14) in ref. [12]Q but 
disagree with their 0(^^) result (the 0(^^) is not displayed in [12]). Notice that for ^ large 
enormous cancellations occur in formula (3) and hence the analytic expansion (4) may prove 
very useful. 



3 Wave Function at the Origin 

The correction for the wave function at the origin for n = 1 states reads 



where \l/io(x) is the Coulomb wave function. The first bracket corresponds to the zero photon 
exchange and has already been calculated analytically in [5] . The second and third brackets 
correspond to the pole subtraction and multi-photon exchange contributions respectively. 
Fi{^) , i = 2,3 are defined as follows: 



Taking e„ = in that reference and upon correcting an obvious misprint ki k. 




-li^e + ^om) + li^e +e- 2)^3(0 + 1 [^e + y ) log | 



(i-e^Fi(o) + 



(6) 



4 



^2(0 



2 

/ 



d0 log 



"sin 6' + ^" 
sin^ 



de 



1 



sin 6* + ^ 



log 



sin6' + { 
sin^ 



(7) 



F2{i) and ^3(0 can be expressed in terms of Clausen integrals and dilogarithms. We present 
the explicit formulas in Appendix A. The key steps in order to obtain (6) are given in 
Appendix B. 



V, 



vpc 



5(X) 



Fig. 2: Diagram rendering the correction to the wave function at the origin. The double hne is the 
Coulomb propagator of the non-relativistic pair and the star a local (5(x)) potential. 

For ^ large, namely mi» iia, (6) behaves like 



5*io(0) 



^— >CXD 



a 



TT 



*io(0) 



Stt 107 
16^ ^ 225^2 




(8) 



This result must be compatible with the one obtained by integrating out the light fermion 
first and then calculating the electromagnetic potential. In the case m» mi » fia/n (for 
simplicity, we are assuming h — h', m being the mass of the non-relativistic particles) we 
expect that a local non-relativistic effective theory is obtained after integrating out the energy 
scale mi and the associated three momentum scale ^mmi for the non-relativistic particle. 
The leading term in (8) corresponds to the contribution that would be obtained from the local 
term induced by the diagram in Fig. 3. The logarithm in the subleading term corresponds 
to the iteration of two delta function potentials in quantum mechanics (see formula (5.5) in 
[9]). The second delta function is due to the contribution to the electromagnetic potential 
of the dimension six photon operator (see [32]) which arises after integrating out a heavy 
particle 
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Fig. 3: VacTium polarization correction to the decay width at leading order in when ^ — >■ oo. 
Energies and momenta of order nii and yjmmi respectively dominate the graph and hence the 
Coulomb resummation leads to subleading effects. 

For ^ small, namely mi « /la, we obtain 



5*io(o)(e 



a 



0) ^ ^,o(o) 

TT 



5-^-iog|-5e=+o(a 



(9) 



We have made the following checks. For ^ large and small, the leading term of (8) and (9) 
agree with formulas (22) and (23) of ref. [4] respectively (the next-to-leading terms are not 
displayed in [4]). For ^ small we can also compare with known results for massless quarks in 
QCD. We agree with the 0{as) correction of formula (69) of ref. [35]. We have also checked 
that the formula (6) reproduces the numerical results obtained for di-muonium and pionium 
in refs. [11] and [5] respectively, and we also agree numerically with the analytical result in 
terms of a non-trivial integral of ref. [4] . 



4 Applications 



4-1 Exotic Atoms 



We have listed in Table I and Table II the corrections to some energy splittings and to 
the wave function at the origin respectively of simple exotic atoms of current interest. This 
purely electromagnetic corrections must be conveniently taken into account if one wants to 
obtain precise information of the strong scattering lengths from hadronic atoms. 



TABLE I 





5E2i-5Eio 
a{E2~Ei) 


SE'si-SEia 
a{E3-Ei) 


5E2i~SE2Q 
aE2 


pK- 


.44593 


.38629 


-.10453 


pTr~ 


.18103 


.15388 


-.056337 


Pl^- 


.13616 


.11548 


-.044443 
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Vacuum polarization induced energy splittings for some exotic atoms. 

TABLE II 





t _ me 


<5^ph*(0) 


<5ps*(0) 

V"/ 




<5*(0) 


K-p 


.21648 


.34290 


.15454 


.09650 


.59394 


K+K~ 


.28369 


.29837 


.12958 


.08785 


.51581 




.57635 


.19613 


.07285 


.06166 


.33064 




.64357 


.18237 


.06549 


.05741 


.30527 




.73738 


.16627 


.05703 


.05222 


.27552 


7r"'"7r~ 


1.00344 


.13338 


.04052 


.04099 


.21490 




1.32550 


.10793 


.02876 


.03184 


.16853 



Vacuum polarization correction to the ground state wave function at the origin of some exotic 

atoms. 



4.2 T(ls) andtt 



The current calculations of heavy quarks near threshold assume that the remaining lighter 

quarks are massless. This approximation is difficult to justify a priori at least in two cases. 
For the T(ls) system the typical relative momentum mbas/2 ~ 1.3 GeV. [34] is of the same 
order as the charm mass rUc ~ 1.5 GeV. The effects of a finite charm mass in the binding 
energy have been recently quantified in [22] . We give in Table III the size of these effects in 
the wave function at the origin. For the ti production near threshold at a relative momentum 
mtQ^s/'^ ~ 18 GcV. the effects of a finite bottom mass ~ 5 GeV. should be noticeable. 
In order to estimate them, we also show in Table III the size of this effect, both for bottom 
and charm, in the wave function at the origin for the would-be-toponium (Is) state. 

TABLE III 



hh 
it 
it 



^ = 



1.4 



.10 



^ = 



.28 



.088 
.011 
.0019 



Vacuum polarization correction to wave function at the origin in quarkonia. MS has been used. 
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If the corrections were organized in a series of cc^/tt multiplied by numbers of order 1, one may 
conclude that the leading effects of a finite quark mass are: (i) in the T(ls) system for charm 
more important than the next to leading corrections [35]; (ii) in the tt system near threshold 
for bottom (charm) as important as (less important than) the next to leading corrections 
[24] . However, the relativistic corrections do not have the tt suppression and some radiative 
corrections are enhanced by factors of /3o- In practise the next to leading corrections are 
comparable to the leading ones, even for the tt system (see [36] for a discussion). This makes 
the actual size of the finite mass effects smaller than the next to leading order corrections in 
all the cases above. 
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Appendix A 

F2{x) in (7) can be expressed in terms of Clausen integrals. We get 



F2{x) = 



2C/2(arcsin x) — |CZ2(2arcsin x) 



2E 



(-1)* 



k=0 (2fc+l)2 



1.831932 



if X < 1 



if X = 1 



—i Li2{—x — \/x'^ — 1) + i Li2{i (^—x + y/x'^ — l)) — 
-iLi2{-x + y/x^ - 1) + iLi2{-i (x + \/x'^ - l))+ 
+1 (iTT + 4 log(2a:) - 4 log(l + ix- i^Jx^ - 1) - ^\o^{l ^ ix ^ i^J x^ - 1)) + 
+2(f-ilog(^ 



(10) 



2 



^) {^—2i arctan( 



) + 



+ log(l + x — \J x"^ — 1) — log(l + i X — i \J x"^ — \] 



+ log(l + i a; + i sjx^ — X) — log(l + x + \/x^ — \ 



if X > 1. 



Recall that the Clausen integral is defined as 

Cl2{x) -.^^ J d9 log (2 sin -\ ^i^-jx^ -X log(ie-*t) - i Li2(e^^) (11) 

^ ^ 



8 



F^{x) in (7) can be expressed in terms of dilogarithms. We get 



2b 



2b 



2b 



2b 



Xl + a + b)^ 



(l + a + b)^ 
( 2b 



2b 



^. ( 2b \ ^. f (a + b) 

+Ll2 7 7- \+Ll2' ^ ' 



a + b)J \{l + a + b)^ 

where a :— x~^ and b :— ^^11=^. 



Li-) 



2b 

{a-b) 



(a + b) 



+ 



(l + a-b) 



+ 



if X <1, 



F3(l)=2-log2 
F3{x) ^ 



iy/x^ — 1 
-Li2 



2ib 



2ib 



2ib 



2ib 



(1 + a + ib) 



(1 + a + ib) 
( 2ib 



+Lt2 7 — + Lt2 



(a + ib)^ 



+2 arctan 



2ib 
{a + ib) 
(1 + a + ib) 

71 



Liq 



2ib , 
(a — ib) 



(a + ib)^ 



+ 



+ 



if x>l, (12) 



where a :— x ^ and b :— ^'^^ 



In order to make contact with the expressions found in the hterature for the massless hmit 
of (3) the following formula is useful 



V'(2n) - 



{n + l)\(n -l-l)\ (2(n - I - 1 - k) - l)\{2{k + /) + 1)! 



(2n- 1)! 



E 

A;=0 



{n - I - k - 1)P{21 + 1 + k)\k\ 



iP{n + I + 1){13) 



Appendix B 

We sketch here the main steps which lead to our analytic formulas. For the energy shift we 
have to calculate {v — y/l — x^) 



5E„i =< nl\Vypc\nl > 



2aEr, 



{n + l)\ 



'-2n-2l-2 



Sn (n - I - 1)\{21 + 1)\ 



1 

/ 



dx- 



X 



21+1 



x + 



2n 



Vl-x^ (2 + X^) 
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X 



F(-(n-/-l),-(n-/-l),2/ + 2;-). 



(14) 



Since I < n the hypergeometric function above reduces to a polynomial 



F(-(n - ^ - 1), -(n - / - 1), 2^ + 2; .) ^ \ f \ t , ' ^[ , , 

j=o (^ - ^ - 1 - j)! (j + 2/ + l)!j! 



Hence 



(15) 



5E, 



nl 



2aEn "g^ I n-l-1 



j=Q 



3 J 



n + l 



\n-l-j- 



2n-2l-2j-2 



I- 21+23+1 

X / dx —Vl^ (2 + X^) . 

and upon writing 



2n 



1 



X + 



^^2n^l ^ + a. 



(16) 



(17) 



and making the change x — >■ sin^ we obtain (3). 

For the wave function at the origin we have to calculate: 

1 InO >< nOl 



(5^„o(0)^„o(0) = lim < nO|5(x) 



E-^En 



E-H E-Er, 



Vypc\nQ > 



(18) 



Upon using the following representation for the Coulomb propagator [37] 



-1 oo I 

l=Q m=—l 



Gi{x,y,E) = -Aki,{2kx)\2kyye-'^^^y^ £ ^-'i(2M^n't\(2Mr(nO 



^,{n' + l-t!^)V{n' + 2l + l) 



(19) 
(20) 



where A;^ = — 2/i-B, (18) can be split into three pieces 

5^'„o(0) = 5p.^'no(0) + 5.pft^'no(0) + 5mp/^^no(0) 



(21) 



The first piece (pole subtraction) corresponds to the term n' = n in the sum (20) and it can 
be calculated using the formulas given above. The remaining pieces read 
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^ i + V 1 - 



n 



^ \ n'-l / 1 2 \ 

^ ^ F -(n-l),-(n'-l);2;^ .(22) 



Again the hypergeometric function above reduces to a polynomial. For n = 1 it reduces in 
fact to 1 and the sum over n' can be carried out explicitly. We obtain: 



1 2 



a r V [I T \ t 

S.ph'^ioiO) + (5„p,vi/io(o) = -M/io(0) / dv—-^^=^- -j=== + log 



where the first term corresponds to the zero photon exchange and the second one to the 
multiphoton exchange. Again the change of variable v cos 6 and a number of manipulations 
allow us to obtain (6) from the above. 
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